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An effective spin-orbit coupling can be generated in cold atom system by engineering atom-light 
interactions. In this letter we study spin- 1/2 and spin-1 Bose-Einstein condensates with Rashba 
spin-orbit coupling, and find that the condensate wave function will develop non-trivial structures. 
From numerical simulation we have identified two different phases. In one phase the ground state 
is a single plane wave, and often we find the system splits into domains and an array of vortices 
plays the role as domain wall. In this phase, time-reversal symmetry is broken. In the other phase 
the condensate wave function is a standing wave and it forms spin stripe. The transition between 
them is driven by interactions between bosons. We also provide an analytical understanding of these 
results and determines the transition point between the two phases. 



Interaction between matters field and gauge field is of 
central importance in quantum physics. Although atoms 
are neutral and do not possess gauge coupling to real elec- 
tromagnetic field, it has been proposed that a synthetic 
gauge field originated from Berry phase effect can be cou- 
pled to atoms by engineering the interactions between 
atoms and a spatially varying laser field [T]. Very re- 
cently, this scheme has been successfully implemented in 
a Rb 87 Bose-Einstein condensate (BEC) [2]. A straight- 
forward generalization of this scheme can also create a 
non-abelian gauge field to cold atoms [3-5]. This opens 
up a new avenue in cold atom physics, that is to study 
how a coherent matter wave responses to external gauge 
fields, in particular, the non-abelian gauge fields. 

There are already various proposals to achieve a non- 
abelian gauge field. For instance, one can start with a TV- 
pod system where TV atomic internal states are coupled to 
one common state by N different laser fields. It results in 
TV— 1 dark states which play a role as internal pseudo-spin 
degree of freedom [3H6]. Tripod (N = 3) and tetrapod 
(N = 4) setup correspond to spin- 1/2 and spin-1 case, 
respectively. Among all possibly achieved configurations 
of non-abelian gauge fields A, one of the most simplest 
case with A x = a x and A y = a y is equivalent to a Rashba 
type spin-orbit (SO) coupling [H[5]. Very recently, BEC 
with SO coupling has been first realized by NIST group 
using a similar scheme described in Ref. [2] [6]. 

During last a few years, it has been extensively stud- 
ied that the SO coupling in an electronic system can lead 
to a novel state of matter of topological insulator which 
has many exotic physical properties [7 j. Nevertheless, so 
far only few papers have studied the SO effect in a bo- 
son condensate [9j [10] . This letter is devoted to study 
the properties of spin- 1/2 and spin-1 condensate in pres- 
ence of a pure Rashba SO coupling. The model under 
consideration is H = Hq + H lnt , 
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where \P = (^,^) for spin- 1/2 case and \£ = 
(^1,^0,^-1) for spin-1 case. Here we consider a 



quasi- two-dimensional situation and in-plane SO cou- 
pling where k = {k x ,k y } and S = {cr x ,a y }, where a XiV 
are spin-1 /2 or spin-1 representation of Pauli matrices. 
For the interaction term, we note that different pseudo- 
spin components are in fact superposition of atomic hy- 
perfine states, and therefore the interactions between 
them could have a quite complicated form. However, 
as an initial effort to understand this rich system, in 
this work we try to simplify the situation by considering 
a simplify interaction form borrowed from conventional 
spinor BEC, namely, for spin- 1/2 
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we shall also focus on the case with g ± = g 2 > and then 
the interaction can also be rewritten as 
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where ft = + S z = — n^, and Co = gi + #12 an d 
C2 = 9i — 9i2- The cq term is SU(2) spin rotational in- 
variant while the c 2 term breaks spin rotation symmetry. 
For spin-1 case, we shall consider the standard interac- 
tion form [8 
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where n = hi + no + n_i and S = ^ a a a ^^. In both 
cases, Co > and 7 = C2/C0 can be either positive or 
negative. 

The main results of this letter include: (I) A boson 
condensation in k = single particle state is always un- 
stable in presence of SO coupling. (II) The ground state 
has two possible phases: one is named as "plane wave 
phase" (PW) where the ground state is found to be a 
single plane wave; the other is named as "standing wave 
phase" (SW) where the spatial wave function of each spin 
component forms an oscillating standing wave, as pre- 
sented in Fig. [2] and [3j (III) The transition between 
PW phase and SW phase depends on the interactions 
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between bosons. For spin- 1/2 case, it is PW phase if 
7 > and SW phase otherwise; while for spin-1 case, it 
is PW phase if 7 < and SW phase otherwise. (IV) 
We often find long-lived metastable state in which the 
system splits into two domains. In the parameter regime 
of PW phase, it is locally a plane wave state with oppo- 
site wave vector in each domain, and an equally-spaced 
array of vortices plays the role as domain wall. We note, 
unlike BEC in a synthetic magnetic field [2 , the Hamil- 
tonian for non-abelian gauge field considered here pre- 
serves both time-reversal and translation symmetry, how- 
ever, the "PW" phase spontaneously breaks time-reversal 
symmetry, which is unconventional in a bosonic system. 
Moreover, both two phases spontaneously break space- 
spin rotation symmetry. This behavior is fundamentally 
different from the effect of SO coupling in electronic sys- 
tems. We have also verified that other terms not consid- 



ered here, such as linear and quadratic Zeeman field, the 
difference between gi and #2, will only shift the phase 
boundary between these two phases, and will not affect 
the essence of the two phases, as long as their strength is 
relatively not too strong. 

Instability of a condensation on k = state: Without 
SO coupling, the ground state of bosons is a condensa- 
tion in k = state. Hence, it arises the question whether 
the SO coupling k • a will have any significant effect on 
a BEC, in particular, if the SO coupling is weak. Here 
let us take spin- 1/2 case as an example, and first ex- 
pand the boson field operator ^J-(r) = e^b^/y/V 
(V is the volumn). If we still assume a boson conden- 
sation in k = state, i.e. (&<j,k=o) = 4>a , the Bogoli- 
ubov Hamiltonian can be written as H^g — B^H^B^ 
where B^ = (&J k , &t-k,&J k , ^-k) and 
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where = k 2 /(2m) + 2gi\cj) cr \ 2 . We then introduce 
a generalized Bogoliubov transformation B^ = Q^B^, 
where B^ = (/?j k , /3^-k? /^-k) where ft^ is a 4x4 ma- 
trix. To satisfy the commutation relation [/3<rk? /?£/ k /] = 
cW'^kk' and [/3crk, Pa'k'} — 0, ^k has to satisfy the rela- 
tion Q^AQk = A where A is a 4 x 4 matrix ( a * ^ V 
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Suppose i^k can be diagonalized by a Vt^ and the diag- 
onal values of Vl\.H\SI\<l are denoted by A^. A^ satisfies the 
equation Det[i^k — \A] = 0, which gives two solutions 
denoted by A^k and A^k- The Bogoliubov Hamiltonian 
becomes = J^ k a A a k/^. k /^-k- We find that one ex- 
citation branch always has a positive imaginary part in 
a large regime in momentum space, even for infinites- 
imal small SO coupling. Therefore, suppose a BEC is 




FIG. 1: (Color online) An example of real (a) and imaginary 
(b) part of unstable branch of excitation spectrum with SO 
coupling. 



initially prepared in the k = state, as SO coupling is 
turned on, some modes will exponentially grow. There- 
fore, a conventional BEC on k = mode is unstable with 
SO coupling, which arises the question that what is the 
actual ground state of a spinor BEC with SO coupling. 

In addition, we note that the spectrum shown in Fig. [I] 
breaks spatial rotation symmetry. This is also an effect of 
SO coupling. Since the Hamiltonian of Eq. [IJonly possess 
a symmetry of simultaneous rotation of both spin and 
space, namely k x + ik y — >• e %e (k x + ik y ) and a x ± ia y —> 
e±ie (°~x =t i°~y)- However, the Bose condensation locks 
the relative phase between different spin components and 
therefore breaks spin rotation symmetry. Consequently, 
this symmetry breaking manifests itself in real space. 

Numerical Simulation for Spin-1/ 2 Case: We imple- 
ment the mean-field approximation and numerically look 
for the condensate wave function that can minimize 
Gross-Pitaevskii energy using imaginary time evolution 
method. For spin- 1/2 case the Gross-Pitaevskii energy is 
written as 

+ |(I^| 2 + KI 2 ) 2 + |(I^| 2 -KI 2 ) 2 } (5) 
What we find are shown in Fig. [2] For 7 = c^jc^ > 
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0, the densities of both component have no particular 
structure, while the phase of both components behave 
as a plane wave, as shown in Fig. [5|a) and named as 
"PW" phase. Time reversal symmetry is broken in this 
phase. For 7 < 0, the wave functions of both components 
behave as standing waves and lead to periodic density 
modulation of both components, as shown in Fig. [2^b) 
and named as "SW" phase. Time reversal symmetry is 
preserved. Moreover, the higher density regime of spin- 
up component coincides with the lower density regime 
of spin-down component, which represents a microscopic 
phase separation, and also represents a spin stripe state. 

Though in the numerical simulation, we have included 
a very weak harmonic trap which helps to avoid artifact 
from a sharp boundary and also simulates the practi- 
cal situation in cold atom experiment, the results can 
be understood from a homogeneous case. With SO 
coupling, the single-particle spectrum becomes E±^ = 
(ft 2 k 2 ±2ttft|k|)/(2ra), where ± denotes different helicity 
(spin parallel or anti-parallel to wave vector). The single 
particle ground state is in the negative helicity branch 
with |k| = K/h, and the wave function is given by 



(6) 



where cp^ = arg(/c a; +ik y ) and </?_k = (pk + n. Let us first 
consider a simple case that assumes the condensate wave 
function is a superposition of two opposite wave vector 
state as 
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Using this ansatz to minimize the interaction energy, it 
is easy to find, for C2 < 0, it favors ol\ = ct2 — 1/V% 
therefore, <p^ ~ coskr and <p^ ~ isinkr. While for C2 > 




FIG. 2: (Color online) Numerical results for spin-1/2 case, 
al and a2 show the phase of condensate wave function of 
both spin-up (al) and down (a2) component increase from 
— 7r (grey regime) to ty (dark regime) periodically in the PW 
regime; bl and b2 show the density of both spin-up (bl) and 
down (b2) component oscillate periodically in the SW regime. 
7 = C2/C0 = 0.4 for al and a2, and 7 = 0.1 for bl and b2. 



0, it favors the case gl\ — 1, ct2 = or ol\ = 0, 0L2 — 1 5 
namely, the wave function is a single plane wave. 

The next question is whether there will be more than 
one single k state, or a pair of {k, — k} states entering the 
condensate wave function ip. In general, one shall assume 
a superposition of all states in the degenerate circle 
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where the amplitude of k is fixed at k/H to minimize the 
single particle energy. For instance, if is independent 
of the angle of k, one can obtain 
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where 9 is the angle of r. This is a symmetric skyrmion 
solution, which has also been proposed by Ref. [TP] . 
However, if one substitutes the ansatz Eq. [7] into the 
energy function Eq. [5] and minimize energy with respect 
to all c^k, we can find the most favorable solution is al- 
ways that ak is non-zero either for a single k or for a pair 
of {k, — k}, and we do not find a parameter regime in 
which the condensate wave function contains more than 
two wave vector components. 

Numerical Simulation for Spin-1 Case: We now move 
to study spin-1 case, whose energy functional is given by 
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The results from numerical simulation are displayed in 
Fig. [3] As shown in Fig. [3j for 7 > it is "SW" phase 
and for 7 < it is "PW" phase. 

The result of spin-1 can be understood with similar 
analysis above. The wave function for single particle 
ground state is now 
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We shall also consider a superposition as ai^k + <^2</>-k- 
Co term is independent of ol\^ however, for ol\ = c\2 = 
l/y/2, (S) 2 = and it is favored when C2 > 0. This state 
is also a special case of so-called "polar" or "nematic" 
phase in the discussion of spin-1 BEC [8]; while for a\ = 
0, ct2 = 1 or ai = 1, ot2 = 0, (S) 2 = 1 and it is favored 
when C2 < 0. This is also called "ferromagnetic" phase. 

Domains and domain wall: Our numerical simulation 
also finds long-lived met ast able states with domains. For 
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FIG. 3: (Color online) Numerical results for spin-1 case, al- 
a3 are density of 1, 0, and —1 component in "SW " regime; 
and bl-b3 are phase of and b3 are 1,0, and —1 component in 
"PW" regime. 7 = 02/ Co = 0.2 for al-a3, and 7 = —0.2 for 
bl-b3. 



instance, if we start with a random initial configuration, 
most cases the imaginary time evolution leads to a state 
in which the system often splits into two domains. In the 
"PW" phase, the system is locally a single plane wave 
state in each domain, and the wave vector is opposite be- 
tween two domains. This is very similar to the situation 
of ferromagnetism, where one always finds a ferromag- 
netic state made up with locally magnetized domains. 

In presence of domains, we find that both up and down 
components of the condensate wave function contain an 
array of vortices, as shown in Fig. [4^a). We have checked 
that the vorticity of all vortices are the same, and the 
vortices in different components locate alternately. The 
vortex array plays the role as a domain wall. Consider an 
array of vortices with same vorticity, located at x = nl 
and y = 0, where n are integers, in a uniform superfluid 
the gradient of superfluid phase d0 at (x, y) is given by 
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00 as a function of (x,y) is shown in Fig. [i] (c-d). As 
one can see, when \y\ > Z, d x — >• tt/1 is a constant and 
d y — > 0. Hence the system is locally a plane wave state. 
To minimize the single particle energy, one requires I = 
ttH/k, namely, the vortex line density increases as the 
increase of SO coupling. 

Hence, we have established the conclusion (I- IV) sum- 
marized above within a mean-field theory. In future stud- 
ies we will include quantum fluctuations. Due to the de- 
generacy of single particle ground states, quantum fluc- 
tuation, in particular, the fluctuation of rotation mode, 
may lead to fragmentation. However, it is known that 
fragmented state is usually very fragile for a realistic 
system with large number of bosons and is not stable 
against external perturbations. In this case, effects such 
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FIG. 4: (Color online) (a) Density of up (or down) compo- 
nent with two domains in the "PW" phase, (b) A schematic 
of vortex array as a domain wall between two regimes with 
opposite wave vector, (c-d) Consider array of vortices sepa- 
rated by distance / located in x-axes, d x (c) and d y (d) (in 
unit of l/l) at position (x,y) as a function of x/l is shown 
for \y\ = I (red dashed line), \y\ = 1.5/ (blue dotted line) and 
I y I = 21 (black solid line). is the phase of condensate wave 
function. 



as anisotropy of trapping potential will break the spatial 
rotational symmetry and pin the direction of plane wave 
or density stripe, and prevent the fluctuation of rotation 
mode from restoring the symmetry. Mean-field results 
become more stable and the predictions of this work can 
be verified experimentally very soon. 
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